Crossover between ionic/covalent and pure ionic bonding in magnesium oxyde clusters 
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An empirical potential with fluctuating charges is proposed for modelling (MgO)„ clusters in both 
the molecular (small n) and bulk (n — > oo) regimes. Vectorial polarization forces are explicitely taken 
into account in the self-consistent determination of the charges. Our model predicts cuboid cluster 
structures, in agreement with previous experimental and theoretical results. The effective charge 
transferred between magnesium and oxygen smoothly increases from 1 to 2, with an estimated 
crossover size above 300 MgO molecules. 
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Crystalline magnesium oxyde is a purely ionic com- 
pound in which the Mg z+ and O z ~~ ions carry a charge 
Z around ±2i In gas phase, the oxyde anion 2 ~ is un- 
stable and spontaneously decays into 0~+c due to the 
strong electron-electron repulsion. As a result, the effec- 
tive atomic charge in the MgO molecule is much smaller 
than 2: independent ab initio calculations by Ziemann 
and Castleman 2 and by Recio et al~ found Z ~ 0.8. In 
the intermediate sizes regime, (MgO)„ clusters are thus 
expected to show intriguing properties due to a partially 
covalent character of the chemical bonding. Beyond con- 
densed matter or molecular physics, these clusters re- 
ceived some special attention in the astrophysics com- 
munity, where they have been involved in the nucleation 
process of dust in circumstellar shells around M-stars. 

Despite the vast amount of experimental 2 *^ and 
theoreticalSiiSiLSiSiifliiLiS investigations on neutral or 
charged clusters, the way and the rate at which chemical 
bonding evolves from ionic/covalent at small sizes toward 
purely ionic in the bulk remains essentially unexplored. 
Because the electrostatic field created by the ions does 
not vanish in finite systems, the highly polarizable oxyde 
anion has a rather floppy and deformable outer electron 
cloud, which could be responsible for a partial screening 
of the repulsion between cations. However, the situation 
is complicated by the possible coordination-dependence 
of the charge transferred locally. 

Theoretical studies of (MgO)„ clusters can be effec- 
tively separated into two groups. Ab initio or density- 
functional theory (DFT) based calculations have been 
performed on specific geometries, in a rather limited 
size rangei 2 * 3 *^*** 8 ^ These works predict that small clus- 
ters exhibit cuboid-like shapes similar to NaCl rock- 
salt clusters. The apparent charge transferred, as es- 
timated from Mulliken populations, is indeed size- and 
coordination-dependent, and lies between 1 and 1.5 for 
2 < n < 13f3a£ More empirical methods have also been 
used to predict optimal structures ,2*12*11 Ziemann and 
Castleman^ and, more recently, Roberts and Johnston, 11 
have used the rigid ion model (Born-Mayer + Coulomb 
interactions) potential with two possible values of the 
charge transferred. When a charge Z = 1 is taken, 
cuboids are preferentially found as the most stable ge- 



ometries. For Z — 2, as in the bulk, small clusters 
show instead hollow, fullerene-like structures. The ef- 
fects of polarization have been studied by Kohler and 
coworkers^ using the polarizable ion model due to 
Rittner. 13 Wilson^ investigated MgO "nanotubes" made 
of hexagonal (MgO)3 rings stacked. For this he developed 
a more sophisticated compressible-ion modetii with ex- 
plicit coordination-dependent polarizabilities. 

None of these empirical potentials account for the dif- 
ferent charges transferred in MgO clusters. Only in 
Ref.^]the authors explicitely employed a size-dependent 
value of the charge Z, using an arbitrary law Z(n) — 
(2(n + l)/(£n + 1). C was taken such that the crossover 
n* between ionic/covalent and purely ionic, for which Z 
equals 1.5, occurs approximately at n* = 20. However, 
the low energy structures found by Kohler and cowork- 
ers significantly deviate from regular cuboids or stacked 
hexagons in this size ranged An improved treatment of 
electrostatics and charge transfer is provided by fluctuat- 
ing charges (fluc-g) potentials**^*!! based on the prin- 
ciple of electronegativity equalization. 18 Such potentials 
have been used in simulations of water^*^* 2 *! and molten 
salts* 2 ! and have recently proven valuable in describ- 
ing the heterogeneous bonding in coated fullcrcncs and 
nanotubes ■ 2 - 2 ' 2 - They have been extended to include dipo- 
lar terms and the corresponding polarization* 2 *^ In MgO 
clusters atomic polarization cannot be neglected, and we 
provide here a self-consistent treatment of these effects. 
Briefly, the system is made of N magnesium cations and 
M oxygen anions, each carrying a charge qi and located 
at the position vector r^. The total potential energy V 
of the system is written as V = V Tcp + Vq. The repulsion 
interaction V rcp is taken in the usual Born-Mayer format 
as a function of the distance between ions i and j: 
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The total electrostatic energy Vq is expressed as 
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£i = £Mg or £o are the electronegativities associated with 
each atomic element. Jy is the Coulomb integral between 
ions i and j, taken in the Ohno representation as^ 
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OLi = «Mg or ao are the atomic polarizabilities, and E, 
is the electric field vector felt by ion i: 
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Finally, the last term in Eq. includes a Lagrange mul- 
tiplier A, which accounts for the conservation of the total 
charge Q of the system. Given an instantaneous set of 
positions {r^}, the charges {qi} are found by minimizing 
Eq. J2J above. Thanks to the linear dependence of the 
electric fields Ej's on the charges, the expression of Vq is 
quadratic in the g,'s, and its minimization can be done 
readily using linear algebra. The charges are solution of 
the matrix equation CX = D, where X = {qi,X} is a 
N + M + 1 vector and D has components Di = —Ei for 
i < N + M, and Djv+m+i = Q- The element of 
matrix C is expressed as 
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with J a — f/R. The polarizability of the cluster can 
be calculated by imposing external electric fields and by 
computing the variations of the electric dipole. The pre- 
vious expressions need to be modified accordingly to in- 
corporate this effect. This model has 11 independent 
parameters, including D, (3, the a^'s, the U^'s and 7y's. 
Only the difference in electronegativities Ae = eug — £o 
is physically relevant; 17 In order to make the model trans- 
ferable from the molecular range up to the bulk, some 
constraints must be imposed on these parameters. In the 
MgO diatomics, the equilibrium distance, charge trans- 
ferred and electric dipole are known (see for instance 
Ref. ^3 an d references therein). For this molecule, one 
must minimize the energy function V with respect to 
the Mg-0 distance r. After some calculation we find 
V M &°{r) = De-P r +VQ g °{r) with the electrostatic term: 



Q 



1 



(Ae) 



2 (a Mg + «o)[^ g0 M] 2 + 2J Mg oM-l^ gMg -[/° ( o ) 

were we have employed the notation J^go = dJugo/dr. 

In the (Bl) crystal, the charge transferred Z and lat- 
tice constant a are chosen as reference data. Using the 
assumption JMgo(a) ~ I/O) one finds the binding energy 
per ion: 
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with M the Madeling constant. For a set of parameters, 
the total energies corresponding to the diatomics and to 
crystal must be minimized with respect to r or a, respec- 
tively. This is done numerically using the Ohno form of 
the Coulomb integral, Eq. © above. The full parameter- 
ization of the model can then be achieved by minimiza- 
tion of an error function x 2 > to reproduce some dimcr 
and crystal properties. The following values have been 
adopted: D = 6056 eV, j3 = 4.89 A" 1 , U^ Mg = 0.46, 

7oo = 

2 



u oo 



L.13, C/^go = 0.85, 7M gMg = 0.35 A~ 2 , 
0.49 A" 2 , 7M g o = 0.36 A~ 2 , Ae = 0.935, a Mg = 0.18 A 3 , 
and ao = 4.65 A 3 . 
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FIG. 1: Lowest energy structures of (MgO) n clusters from 
Monte Carlo minimization using the self-consistent polariz- 
able fluc-q model. 
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This set predicts the charge transfer to be 0.92 in the 
MgO molecule, and 1.92 in the crystal. The equilibrium 
distance is 1.86 A in the molecule and the crystal lattice 
constant is 1.78 A. These data correspond to respective 
errors of 6% and 16% when compared with the experi- 
mentally measured values. The binding energies cannot 
be compared with reference values, because the present 
model includes extra self-energy U° terms. 

The lowest energy structures of (MgO)„ clusters have 
been searched using the basin-hopping or Monte Carlo 
minimization algorithm^£*2& For each size in the range 
2 < n < 30, 5000 Monte Carlo steps were performed 
starting from a random guess geometry. We also locally 
optimized databases of structures found by global opti- 
mization of the rigid ion model with fixed charges ±1 and 
±2. In many cases, the global minimum was found to lie 
within the database obtained with Z — ±1. 

The structures of the global minima are represented 
in Fig. ^ Beyond n = 3, they are based on small 
(MgO)4 cubic units. The cuboid picture found in the 
present Communication is essentially similar to the re- 
sults of Roberts and Johnston 11 obtained with a genetic 
algorithm, except for the slight distortions due here to 
polarization. We do not find evidences for hollow^iil 
or "spiky" geometries^ Stacking of (MgO)3 hexagonal 
units leads to isomers slightly less stable than cuboids 
at the same size. This partly explains why the global 
minima at n — 14 and n = 22 differ from the results of 
Roberts and Johnston>iI Actually the fact that hexag- 
onal rings are less favored in the present model is not 
in strong contradiction with experimental results, for 
two reasons. Firstly, experiments have been performed 
by mass spectrometry on charged species, which may 
well exhibit different structures than neutrals. Secondly, 
most magic numbers peaks interpreted as the signature 
of hexagonal stacks are indeed compatible with cuboid 
like geometries. 

The variations with size of the binding energy of the 
lowest energy structures found with the present polar- 
izable fluc-g model are depicted in Fig. [5] The bind- 
ing energy E shows a global increase, which can be fit- 
ted approximately in a liquid drop fashion as E(n) — 
-15.13n + 1.83n 2 / 3 + 0.8971 1 / 3 + 1.45. The values of the 
latter parameters are slightly changed if we include larger 
cubic clusters such as (MgO)iog- By construction, the 
crystal binding energy found from this expression is close 
to the numerical minimization of Eq. Q . 

As can be noted in Fig. [3 there are some devia- 
tions from the smooth behavior of the fitted energy. 
To see them more clearly, the second energy difference 
A 2 -E(n) = 2E(n) — E(n + 1) — E(n — 1) has been repre- 
sented in the inset of this figure. This quantity is usually 
convenient to find the special stabilities of some sizes. 
The most stable clusters appear here at n = 2, 4, 6, 9, 
12, 18, and 24. The sizes n = 15, 21, and 27 can be 
added as relatively stable. All these clusters are perfect 
cuboids. The magic character of the n = 15 cluster is less 
marked, due to the fact that (MgO)i6 is also a cuboid. 
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FIG. 2: Binding energy of (MgO)n clusters in the range 
1 < n < 30. The horizontal dashed line marks the asymptotic 
bulk limit. The open circles correspond to the structures of 



Fig-H the solid line is a fit of the form E(n) jn = a + bn 



-1/3. 



cn 2,/3 + dn i , with parameters given in the text. Inset: sec- 



ond energy difference A2 E(n) = 2E(n) 
versus n. 
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We turn now to the problem of charge transfer, and 
more generally to the ionic or covalent nature of the 
chemical bonding in MgO clusters. In Fig. we have 
represented the average charge (q) carried by the ions in 
the cluster, regardless of their position inside the clus- 
ter or coordination number. This quantity is defined 
as the mean value over all magnesium ions, which is 
the exact opposite of the mean value over all oxygen 
ions. This definition is somewhat loose and arbitrary, 
because all ions do not play the same role in the clus- 
ter due to the large surface/ volume ratio. From Fig. [31 
we see that charge transfer is strongly size-dependent in 
the present fluctuating charges model, and that the con- 
vergence toward the bulk limit is much slower than as- 
sumed by Kohler and coworkers. !£• In fact, knowing that 
the average charge effectively reaches about 2 electrons 
at large n allows us to fit the variations of (q) with n 



as (q)(n) 



an 
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The effec- 



tive crossover size between mixed ionic / covalent bonding 
and pure ionic bonding can then be estimated as n* such 
that (q)(n*) ~ 3/2. To get more realistic values for n*, 
we have included the data corresponding to larger cubic 
clusters, namely (MgO)32 and especially the 6x6x6 
cluster (MgO) io8- Because the size range covered re- 
mains relatively small, the value for the latter cluster 
was given a weight of 10 in the fitting process, with re- 
spect to smaller sizes. Using this procedure, we find the 
crossover size to be located at n* ~ 300 ± 100 depending 
on the presence of the large cluster in the fit. Including 
the value for (MgO)ios results in an increase of n* , and 
including the data for larger clusters should further shift 
the crossover size toward several hundreds or thousands 
MgO molecules. 
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Within the present empirical model, electrostatic prop- 
erties are naturally coordination dependent. In the inset 
of Fig.|3| the modulus of the charge transferred is plotted 
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FIG. 3: Average modulus of the charge transferred over 
all ions in (MgO) n clusters. The lines are fits of the form 
(q){n) = 2 - a'n~ 1/3 - b'n~ 2/3 + en' 1 , with (solid line) or 
without (dashed line) the value for the (MgO)ios cubic clus- 
ter. Inset: effective moduli of the charges carried by all mag- 
nesium (full squares) or oxygen (empty circles) ions in the 
(MgO)io8 cluster versus their coordination number. The lines 
are a guide to the eye. 

for each ion in the (MgO)ios nanocrystal versus its coor- 
dination number. Several features are of interest. First, 
the magnitude of the charge transferred increases with 
coordination, as expected from the decreasing intensity 



of the electric held. This is in agreement with the elec- 
tronic structure calculations performed by Recio et al., 3 
by Veliah et al.fi and more recently by Coudray et al~ 
Second, the charge carried by magnesium ions is more 
sensitive to coordination than the charge carried by oxy- 
gen ions. This is also in agreement with the findings of 
Veliah et alfi The above results confirm that conventional 
potentials with fixed charges are not fully appropriate 
to describe MgO clusters. This had been addressed by 
Wilson who considered phenomenological coordination- 
dependent polarizabilities within the compressible-ion 
model?i2, Coordination-dependent charges are a natural 
outcome of the present potential, allowing to study MgO 
clusters in a wide range of condensed phases. While the 
present potential is able to treat large clusters beyond 
the possibilities of first principles-based computations, 
the needed inversion of a square matrix can be a limiting 
factor. Fortunately, extended Lagrangian techniques 1 ' 
can reduce the computational cost significantly, making 
the polarizable fluc-g model valuable in thermodynamical 
context. 

To conclude, we proposed an empirical model to de- 
scribe ionic/covalent bonding in MgO clusters. This 
model is based on fluctuating charges and incorporates 
atomic polarization in a self-consistent way. By fitting 
it on both molecular and bulk properties we found that 
small clusters preferentially exhibit cuboid geometries, 
showing magic numbers in good agreement with experi- 
ments. The average charge carried by magnesium or oxy- 
gen atoms smoothly increases, and the crossover between 
ionic/covalent and pure ionic bonding was estimated to 
be above 300 molecules. The model correctly predicts 
that the charge transferred depends on coordination. 
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